%% Producing Figure 3

close all
warning off

load('no_policy')

%% plot histogram for pr and pd

close all

figure('Units','Inches','OuterPosition',[1 1 7 7])
hist  = histogram(pr_next(1000:end),100,'Normalization','Probability','FaceAlpha',0.5,'LineWidth',1,'EdgeAlpha',0);
%hist  = histogram(pr_next(1000:end),(0:0.005:0.5),'Normalization','Probability','FaceAlpha',0.5,'LineWidth',1,'EdgeAlpha',0);
%ylim([0 0.9]);
%xlim([0 0.5]);
%xticks(0:0.1:0.5)
grid on
set(gca,'FontSize',30);
set(gca,'yticklabel',num2str(get(gca,'ytick')','%4g'))

figure('Units','Inches','OuterPosition',[1 1 7 7])
hist  = histogram(pd_next(1000:end),100,'Normalization','Probability','FaceAlpha',0.5,'LineWidth',1,'EdgeAlpha',0);
%hist  = histogram(pd_next(1000:end),(0.5:0.005:1),'Normalization','Probability','FaceAlpha',0.5,'LineWidth',1,'EdgeAlpha',0);
%ylim([0 0.9]);
%xlim([0.5 1.0]);
%xticks(0.5:0.1:1)
grid on
set(gca,'FontSize',30);
set(gca,'yticklabel',num2str(get(gca,'ytick')','%4g'))
